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ANALYSIS  OF  HBUCOUM.  nUK  AD  CHBQCAL  REACTION 
1*  A  KUUfC  LAYER 


by 


B.M.  Cetegan  and  V.A.  Slrlgnano 
University  of  California 
Irvin* ,  California  92717 


This  paper  concerns  an  analytical  study  of 
molecular  nixing  and  flnlta  rata  chaaical 
raactlona  in  an  Inflnit*  row  of  vorticaa 
representing  a  nixing  layar.  Diffusion  aquations 
for  raactlng  and  non-raactlng  spaclas  ar*  solvsd 
locally  In  tha  Lagranglan  frasn  of  rafaranca 
following  notarial  alaaanta.  Tha  flovflald  la  a 
given  for  tha  problaa .  Tha  concantration 
distributions  In  tha  vortax  structurs  ar*  composed 
froa  thaaa  analytical  solution*  and  presented  for 
several  cases.  Tha  probability  density  functions 
constructed  froa  those  distributions  are  coaparad 
with  tha  alxlng  layar  experlaanta.  For  tha  cases 
studied,  the  coaq>arlsons  Indicate  encouraging 
agraaaant  for  tha  raactlng  flow  casa.  However, 
soaa  dlffarancas  ar*  present  for  the  non-raactlng 
flows. 


In  aany  coabustlon  systaaa,  tha  combustion  of 
fual  and  oxidizer  takes  place  in  turbulent  flew 
reactors  where  fuel  and  oxldlser.  In  either 
prealxed  or  non-prealxad  form,  ar*  Introduced  into 
a  alxlng  sons.  Sine*  chealcal  reactions  taka 
place  upon  aolecular  alxlng  of  reactants,  there  la 
a  strong  Interaction  between  tha  coabustlon 
processes  and  the  local  flowflsld.  A  typical 
oxaapl*  la  tha  reacting  alxlng  layar  behind  a 
splitter  plat*  separating  tha  fual  and  oxldlcar 
streams  that  aeet  at  tha  plate' a  trailing  edge. 
In  general,  the  two  straaaa  have  different  speeds 
to  preaot*  alxlng  of  the  reactants ,  thus  tha 
coabustlon  rat*.  In  such  a  flow  configuration, 
well-known  Initially  two  dlaanalonal  vortices  are 
formed  due  to  the  onset  of  tha  Kelvln-Helaholtz 
Instability.  This  problaa  ha*  racalvad  vast 
amount  of  attention  In  tha  past  several  dacadas. 
Experimental  studies  on  non-raactlng  and  raactlng 
shear  layers'1'*'3'*'  have  bean  complemented  by 
elaborate  numerical  simulations  of  this  flow 
configuration'3'*'''.  The  reference*  cited  hare 
ar*  Just  a  few  exaaplas  of  aany  in  tha 
literature.  These  studies  have  made  significant 
contributions  to  our  unde -standing  of  these 
technologically  Important  flowflalds. 

In  maerlcal  simulation  of  reacting  flows, 
computations  aust  be  performed  that  resolve  saall 
length  scales  (of  the  order  of  the  diffusive 
length  scales)  to  predict  nolacular  alxlng  and 
chealcal  reactions.  This  level  of  resolution,  In 
general,  requires  very  time  consuming  and 
expensive  calculation*.  At  tha  same  tin*,  some 
aathode ,  such  as  finite  difference  simulations , 
ar*  hindered  by  the  numerical  diffusion  errors 
which  need  to  b*  eliminated  by  elaborate  and 
computer  time  Intensive  numerical  methods. 

The  problaa  considered  bar*  1*  an  approximate 
analytical  approach  to  solve  tha  alxlng  problaa  In 


a  alxlng  layer  labadded  with  a  line  of  two 
dimensional  vortices.  This  situation  approximates 
the  taaporal  growth  of  a  alxlng  layer  locally 
without  the  effect  of  pairing.  In  this  paper,  the 
alxlng  and  reaction  of  two  initially  segregated 
species  streams  are  analyzed  In  a  Lagranglan  frane 
of  reference  accounting  for  the  stretching  of  the 
notarial  lines  as  aolecular  diffusion  and  chemical 
reactions  are  taking  place  across  them.  This 
approach  follows  the  methodology  of  Marble'8'  and 
Karbl*  and  Karagozlan'*'  In  that  the  study  of  the 
local  behavior  does  provide  an  accurate 
representation  of  the  local  physical  processes. 
This  approach  was  recently  used  by  Cetegen  and 
Slrlgnano'10'  to  determine  the  concentration  field 
of  non- reacting  and  reacting  scalars  vlthln  a  two- 
dimensional  single  vortex  and  Its  mixing 
atatistics.  In  this  paper,  this  technique  Is 
extended  to  the  study  of  an  Infinite  row  of 
vortices  which  exhibit  some  differences  with 
respect  to  the  strain  field. 

In  the  following,  the  aethod  of  analysis  Is 
first  outlined.  The  approximate  treatment  of  the 
finite  rate  chealcal  reactions  between  the  two 
reactive  species  Is  described.  The  concentration 
distributions  and  alxlng  atatistics  are  presented 
and  discussed  In  the  results  section. 


METHOD  OF  ANALYSIS 


Flowflsld 

The  computation  of  the  concentration  field 
within  an  Infinite  row  of  vortices  relies  on  the 
knowledge  of  the  flowflsld.  The  velocity  field  Is 
utilized  to  calculate  tha  successive  positions  of 
the  material  element*  In  tine  as  well  as  the 
strain  rates  experienced  by  them.  This 
facilitates  analytical  solution  for  the 
concentration  profiles  about  the  material  elements 
locally.  In  this  paper,  the  velocity  distribution 
for  an  Infinite  row  of  point  vortices  Is  aodlfled 
to  Include  a  viscous  core  In  the  center  of  each 
vortex.  The  velocity  distribution  for  an  infinite 
row  of  point  vortices'11'  Is  given  by 
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The  coordinate  system  and  the  parameters  in  these  " 
equations  ar*  depicted  In  Figure  1.  The  velocity 
distribution  represented  by  Equation  <1)  -  - 

represents  a  horizontal  velocity  Jump  of  too/ 

magnitude  T/a,  froa  -T/2a  to  +r/2a.  Since  the  ^ , y  - 

tangential  velocity  at  the  center  of  a  potential 
point  vortex  tends  to  Infinity,  a  situation  which  1  anA/or 
Is  not  realised  In  viscous  fluids,  wo  can  *ci£i. 

i  i 


□  □ 


Incorporate  the  affacta  of  viscosity  by 
Introducing  •  tiao  dependent  viscous  cora  given  by 
a  Banklne  vortex. 

tinea  tba  problem  is  periodic  in  tba  x 
direction,  w  can  consider  one  of  the  vortices 
under  the  influence  of  seal -infinite  rows  of 
vortices  on  each  side.  In  particular,  if  ve 
confine  our  attention  to  tbs  vortex  situated  at 
x  —  0,  y  -  0,  then  tba  velocity  field  can  be 
obtained  by  superposition  of  a  viscous  core 
yielding. 
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can  be  invoked.  Here,  c(t)  is  the  tine -dependent 
strain  rate  of  the  aleaent  along  Its 
tangent,  f  and  r  are  the  stretched  transverse 
distance  and  the  transforaed  tine  accounting  for 
the  strain  history  of  that  particular  eleaent 
until  it  arrives  to  Its  current  location.  Under 
this  transformation,  Equation  (3)  reduces  to  the 
ona-dlaenslonsl ,  unsteady,  diffusion  equation 

ac  _  a2c 

a7  -  D  T-2  <5> 

of 

with  the  appropriate  initial  and  boundary 
conditions.  For  the  case_  of  C(y)  >  0,  t-0) 

-  1.0,  C(y  <  0,  t-0)  -  0.0,  C(y  v  »,t)  -  1.0  and 
C(y  t)  -  0.0,  the  solution  can  be  written 

in  terns  of  the  error  function  as 


c  -  f  [1  ♦  erf  (-ir:)] 


Sin  (^) _ 
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The  superposition  is  valid  in  the  potential  flow 
region  outside  the  viscous  core  end  is  approximate 
In  the  core  region.  However,  this  approxiaation 
is  aduissable  in  view  of  the  fact  that  aost  of  the 
straining  occurs  in  the  potential  region  outside 
the  viscous  core. 

Molecular  Diffusion 

In  this  flowfleld,  a  aaterlal  surface 
initially  coincident  with  the  horizontal  x-axls, 
undergoes  a  distortion.  Across  this  aaterlal 
surface  separating  the  two  species,  aoleculsr 
diffusion  takes  place  as  this  interface  is 
stretched  and  transported  in  the  vortex  field. 
The  species  diffusion  equation  for  a  locally 
planar  Interface  can  be  written  as 

dx  dy  9y*  a* 

Here.js  and  y  ere  the  local  orthogonal  coordinates 
with  x  along  _  the  tangent  of  the 
eleaent,  u  sad  v  ere  the  corresponding  velocity 
eeupousKts,  C  is  the  species  aass  fraction,  and  D 
is  the  biliary  diffuslvlty  which  is  assumed  to  be 
egua£  for  both  species.  The  local  flowfleld 
(u,  v)  can  he  directly  related  to  tba  local  strain 
field  experienced  about  the  materiel  surfaces. 
The  strain  rates  con  be  obtained  from  the 
flowfleld  given  in  Equation  (2).  If  this 
flowfleld  produced  a  spatially  uniform,  tlae- 
dspowdeat  normal -strain  field,  than  the  second 
terns  ea  both  sides  of  Equation  (3)  could  be 
eliminated.  Furtbaraore ,  a  aaterlal  transforms- 
tlea  given  by 

-  Jl  «**> 

r  -  ye 

e  -  J*  (exp  <2j*'«(f)  df)  d f  <*) 


In  calculating  the  strain  rates  _in  the  local 
orthogonal  coordinate  system  (x,y)  attached  to 
the  aaterlal  element,  it  is  recognized  that  the 
local  flowfleld  produces  both  a  normal  and  a  shear 
strain.  While  the  normal  strain  component  has  a 
strong  Influence  on  the  molecular  diffusion  across 
the  material  surface  through  steepening  of  the 
concentration  gradients,  the  effect  of  shear 
strain  on  the  diffusion  process  is  much  less 
significant.  Therefore,  the  calculation  of  the 
strain  rates  for  the  solution  of  the  diffusion 
equation  can  be  based  on  the  normal  strain  only. 
Until  this  point,  much  of  the  method  of  analysis 
follows  that  of  Marble'8'  and  Marble  and 
Karagozlan'9' .  In  the  determination  of  the 
concentration  distribution  about  material 
surfaces,  it  is  found  that  the  diffusion  cones 
surrounding  the  neighboring  aaterlal  surfaces 
overlap  in  the  tightly-winded  core  region  of  the 
vortex.  This  effect  becomes  more  pronounced  as 
the  8chaldt  number  (Sc  -  v/D)  decreases  or  the 
aoleculsr  diffuslvlty  Increases.  This  overlap  Is 
taken  into  account  by  superposition  of  the 
concentration  profiles  fro a  neighboring  regions. 
The  superposition  is  adaissable  in  view  of  the 
linearity  of  the  transforaed  diffusion  Equation 
(3).  This  permits  the  construction  of  the  whole 
concentration  field  within  the  vortex  from  local 
concentration  profiles  at  different  locations 
within  the  vortex.  Soae  error  is  introduced  by 
the  spatial  variation  of  the  strain  rate,  t. 

The  concentration  fields  evaluated  in  this 
aatmer  are  subsequently  used  to  determine  the 
probability  density  distributions  (pdf)  fron  the 
concentration  distribution  along  the  x-dlrectlon 
at  different  y-posltlons  utilizing  Taylor's  frozen 
flow  hypothesis. 

The  periodicity  in  the  x-dlrection  allows  the 
construction  of  pdf’s  from  the  concentration  field 
within  a  single  vortex  In  the  row.  In  cosq>arlng 
the  calculated  pdf's  with  the  experiments,  soae 
differetices  are  expected  due  to  the  Lagranglan 
nature  of  the  computations  and  the  Eulerlan  nature 
of  the  experiments. 


Finite  rata  chemical  reactions  within  a 
vortex  are  treated  with  a  aethod  based  on  Creen's 
function  solution  of  the  reacting  species 
conservation  equation.  The  production  terms  in 


eh*  tp*cln  conservation  aquation*  can  b* 
eliminated  by  tha  appropriate  linear  combination 
of  thaaa  aquation* .  Thla  particular  dependant 
variable ,  called  tha  Schvab-Zel'dovich  variable 
glvaa  by 

c.vevv  <7) 

Yol  +  *of  Yfl 

aatlsflae  Equation  (3)  provided  that  the  binary 
diffualvltlea  are  all  equal.  Hera,  Y’s  are  maaa 
fraction*,  4 1*  tha  stoichiometric  oxidant  to 
fuel  Base  ratio,  eubecrlpte,  f,  o  and  1  refer  to 
fuel,  oxidant  and  Initial  raepactivaly.  In  tha 
a  lap la  case  of  a  eacond- order  finite -rata  reaction 
between  Coo  apaclae  (fuel  and  oxidizer),  only  one 
of  tha  species  conservation  aquations  need*  to  be 
solved  and  used  along  with  the  solution  for  the 
Schvab-Zal’dovlch  variable  to  dataralna  completely 
the  reacting  scalar  fields.  Tha  transformed  fuel 
species  equation  under  transformation  in  Equation 
(4)  can  be  written  as 

«Y,  ah.  w  M 

_  D  +  h(0  — i  (8) 

Br  »(2  * 


Here,  h(r)  —  dt/dr,  and 


related 


stretching,  M*  end  p  are  tha  fuel  aolecular  weight 
and  oedluei  density  raepactivaly.  Solution  of 
Equation  (•)  can  be  obtained  by  the  convolution  of 
the  appropriate  Green’s  function'11 ,12) 

2 

C  lt.rir.fi  -ii^e*1**  <,*,-')Coeo  (f  -  f’)do 


axpraselon  with  the  derivative  of  Equation  (11), 
tha  following  two  intagro- differential  equations 
for  f  and  g  are  obtained 


|f  ♦  Do2  f  -  -f  Vf  ’  ’T>  h(r)  Co*  ar'dr' 
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If  +  Da2g  -  J^2  Vf(r.T)  h(r)  Sin  of’dr 


They  are  subject  to  the  initial  conditions 
f(a,0)  -  0  and  g(o,0)  -  0.  In  order  to  iwplenent 
tha  solution,  an  approximation  is  wade  in  Equation 
(11)  by  discretizing  the  continuous  spectrum 
of  a  into  discrete  values  from  zero  to  a 
sufficiently  large  maximum  value  a  .  Equation 
(11)  becomes  " 

rf(t.T)  - 

If  (o  ,r)  Cosa  f(Aa)  til  («  ,r)Sln  o  f(da) 
n  n  n  “  n  n 

n  n 


provided  that  f  and  g„  -»  0  as  o  -»  o  .  Under 
this  approximation.  Equations  (12a)  "and  (12b) 
transform  to  the  following  Integro-dlfferentlal 
aquations : 


Tha  solution  along  a  radial  ray  can  be 
obtained  If  tha  layers  of  reactants  within  the 
vortex  are  treated  as  multiple  strips  of  reaction 
layers  undergoing  straining.  This  permits 
determination  of  the  fuel  mass  fraction 
distribution  along  a  radial  ray 

lf((,r)  - 

Si  S*Z  «f<r'*T,)  Mr’)  C(r,r|f',T’)  df’dr’ 


df  _  H. 

T-®  +  Do  r  -  —  y  VAS'.r)  h(T)  Co*  a  rdf' 
dr  n  n  pn  J  -m  f  '  n’ 


9 

♦  D  "n  V  n  Vf’’°  h(r)  Sln  °nr'df’ 


In  this  expression,  l!  is  the  solution  for  Yf  in 
the  absence  of  chemical  reactions.  Thus,  it 
corresponds  to  tha  pure  diffusion.  Solutions 
obtained  along  radial  lines  are  used  to  construct 
the  field  of  reacting  scalars. 

Solution  of  Equation  (10)  can  be  simplified 
since  the  Fourier  fora  of  the  Green’s  function  in 
Equation  (7)  suggests  a  fora  of  Yf  given  by 

tf(r.o  - 

/I  f(m,r)  Cos  nfda  -  g(o,r)  Sin  of  do 


where  f(a,r)  and  g(e,r)  are  functions  to  be 
determined.  After  differentiating  Equation  (10) 
with  respect  to  r  and  comparing  the  resulting 


with  fn(0)  -  0  and  gn(0)  «  0.  These  equations  are 
Integrated  by  a  second  ordar  Rungs -Kutta  method  to 
determine  tha  amplitude  functions  fn  and  gp. 

The  chemical  reaction  considered  here  Is  a 
second  ordar,  dilute  Isothermal  reaction  between 
HO  and  O3.  The  reaction  products  are  NOj  and 
and  both  species  are  diluted  in  Hj-  Th*  reaction 
rats  is  taken  from  Clvl  st  al'13',  as 


5  .  x 

,  "  -2-2  *  I®’  Vo,  ”C 


The  initial  mass  fractions  of  the  species  were 
taken  to  be  0.01.  In  the  next  section,  the 
results  for  reacting  and  nonreacting  cases  are 
discussed. 
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In  this  section,  computed  spsclss  asss 
fraction  distributions  snd  tbs  corresponding 
probability  density  distributions  within  reacting 
and  non-reacting  row  of  vortices  are  presented  for 
three  cases.  Figure  2  shows  the  distortion  of  an 
initially  horizontal  Interface  separating  two 
species  In  the  upper  and  lower  planes  In  the 
presence  of  an  infinite  row  of  vortices.  In  this 
figure,  the  two  species  are  lsalsclble  (black  and 
white) .  this  vortex  configuration  simulates  a 
shear  layer  between  two  streaas  of  fluids  with 
equal  but  opposite  velocities.  In  our 
coeqputatlons ,  we  consider  a  reacting,  gaseous 
shear  layer  containing  ozone  and  nitric  oxide  as 
reactants.  Both  reactants  are  diluted  In  nitrogen 
so  that  the  chealcal  reaction  la  Isothermal.  For 
the  first  two  cases  studied,  ozone  and  nitric 
oxide  with  a  mass  fraction  of  0.01  each  occupy  the 
upper  and  lower  half-planes,  respectively.  In  the 
third  case,  nitric  oxide  occupies  the  upper  half 
plane  with  a  mass  fraction  of  0.136  and  ozone 
occupies  the  lower  half-plane  with  a  mass  fraction 
of  0.01.  This  case  slaulatas  the  reacting  shear 
layer  experiment  of  Masutanl  and  Bowman*33. 

The  first  case  corresponds  to  s  non- 
dimensional  vortex  strength  (or  vortex  Keynolds 
number)  r/2«v  of  50,.  a  dimensionless  time 
paraamter,  r  —  4vt/a/  of  0.02  and  a  Schmidt 
number,  »/D  of  unity.  Figure  3  shove  the 
normalised  maze  fraction  of  ozone  acroas  the 
vortex.  The  reacted  core  Is  surrounded  by  some 
excess  ozone  In  the  miter  layer  of  the  vortex  as 
seen  in  Figure  2 .  The  corresponding  M0 
distribution  Is  shown  In  Figure  4.  It  can  be  seen 
that  soma  MO  la  still  present  In  the  reacted  core 
because  of  the  non-unity  stoichiometry  of  the 
reaction  (l.e.  stoichiometric  0*  to  MO  mass  ratio 
Is  1.6).  Figure  3  shows  the  distribution  of  the 
conserved  scalar  (mass  fraction  of  the  upper  plane 
species  for  non-reacting  case  or  the  Schvab- 
Zel'dovlcb  variable  for  reacting  case)  In  the 
vortex.  In  the  absence  of  chealcal  reactions,  the 
vortex  assumes  a  more  diffused  structure  with 
lower  concentration  gradients  compared  to  the 
reacting  vortex.  Flguree  6  and  7  compare  the 
probability  density  functions  (pdf)  for  non¬ 
reacting  end  reacting  cases  respectively.  For  the 
non- reacting  vortex,  the  pdf  consists  of  two  peaks 
at  the  unmlxed  concentrations  (Y  -  0.0  and  Y  ■ 
1.0)  as  well  as  finite  probabllltes  In  the  wall- 
mixed  core  region  around  Y  -  0.3.  Chealcal 
reactions  alter  the  pdf  by  shifting  the 
probability  peaks  at  mixed  concentrations  to  lower 
concentrations  and  by  modifying  the  pdf  peak  at  Y 
—  0.0.  Both  of  these  effects  are  due  to  the 
depletion  of  the  species  by  finite  rate  chealcal 
reactions . 

The  second  case  corresponds  to  a  non- 
dimensional  vortex  strength  of  123  and  a  time 
parameter  of  0.004.  These  parameters  simulate  the 
ceodltiona  In  the  early  part  of  reacting  shear 
layer  of  Masutanl  and  Bowman*33 .  The  vortex 
strength  was  computed  based  on  the  reported 
experimental  conditions.  The  time  parameter  was 
estimated  from  the  mean  velocity  at  which  the 
vortices  are  eonvocted  downstream  and  the  location 
at  which  the  comparisons  were  to  be  made  (x  -  11 
cm).  The  ozone  mass  fraction  distribution  Is 
depleted  In  Figure  B.  The  concentration  field  is 
lose  diffused  and  shows  higher  concentration 
gradients  compared  to  the  previous  ease.  This  Is 


appropriate  for  a  vortex  with  higher  clrculetlon 
at  an  early  time  In  Its  growth.  The  higher  vortex 
strength  Increases  the  concentration  gradients 
through  stretching.  Figure  9  shows  the  HO 
concentration  distribution  within  the  vortex.  The 
qualitative  features  of  the  distribution  are 
similar  to  the  previous  case.  The  conserved 
scalar  field,  shown  In  Figure  10,  appears  to  be 
somewhat  different  with  steeper  gradients. 
Therefore,  It  shows  the  vortex  structure  much  more 
clearly. 

The  corresponding  pdf's  for  reacting  end  non- 
reacting  vortices  are  shown  In  Figures  11  and  12 . 
respectively.  Similar  to  the  previous  case,  the 
effect  of  chemical  reactions  on  the  pdf  Is  to 
modify  the  probability  peak  topography  at  zero- 
concentration  of  the  reacting  species  and  to 
reduce  the  probability  peaks  at  the  mixed 
concentrations.  Figures  12  and  13  compare  the 
computed  and  measured  pdf's  for  the  non- reacting 
shear  layer  experiment  of  Masutanl  end  Bowmen 
Although  the  shape  of  the  peaks  at  C  -  0.0  and  1.0 
are  similar,  the  computations  do  not  exhibit  the 
experimentally  observed  ‘hump*  on  the  high  speed 
aide  of  the  mixing  layer.  In  the  Lagranglan  frame 
of  reference  In  which  the  computations  are 
performed,  the  concentration  field  and  the  pdfs 
are  antisymmetric  about  the  x-axls.  Therefore, 
these  calculations  are  not  expected  to  reproduce 
the  ‘hump*  observed  In  the  experiments.  This 
experimental  observation  Is  attributed  to  the 
higher  rate  of  entrainment  of  the  high-speed 
stream  according  to  Dlmotakls*153 .  After 

relaxation  of  the  frozen  flow  approximation,  the 
antlsyi  trie  character  may  no  longer  exist  end 
this  phenomsnon  can  be  studied  theoretically. 

In  the  measured  pdfs  reported  by  Masutanl  and 
Bowman™3  for  the  reacting  mixing  layer,  the  high¬ 
speed  scream  (upper  plane)  and  the  low- speed 
stream  (lower  plane)  contained  NO  and  0j 
respectively.  The  MO  to  Oj  mole  fraction  ratio  In 
the  free  streasu  was  0.11B  for  one  of  the  cases 
they  studied.  Coeducations  were  performed  for 
this  configuration  to  coaq>sre  the  pdf  for  the 
chemically  reacting  mixing  layer.  Figure  14  shows 
the  computed  NO  mass  fraction  In  the  vortex. 
Figures  IS  and  16  give  a  comparison  of  the 
experimentally  determined  and  computed  pdfs  for 
the  reacting  species  0j.  It  should  be  noted  that 
the  shapes  of  the  probability  peaks  are  quite 
similar.  Penetration  of  NO  from  the  upper  plane 
(high- speed  stream)  to  the  lower  plana  (low- speed 
stream)  is  a  result  of  the  abundance  of  NO  In  the 
upper  half-plane.  This  causes  the  extension  of 
the  peek  at  Y  -  0.0  to  the  lower  half -plana.  The 
same  trend  Is  found  in  the  computations.  It  Is 
also  Interesting  to  note  that  the  ‘hump*  observed 
for  the  non- reacting  mixing  layer  has  disappeared 
for  the  reacting  case.  This  Is  believed  to  be 
simply  due  to  the  depletion  of  NO  by  the  chemical 
reaction. 


Computation  of  species  concentration 
distributions  In  non-reacting  and  reacting  row  of 
vortices  have  been  preaented  for  a  few  cases.  It 
was  shown  that  a  piecewise  analytical  approach  to 
the  solution  of  molecular  mixing  and  finite  rate 
chealcal  reactions  can  be  implemented  to  solve  for 
the  concentration  fields.  This  method,  first  used 
by  Marble**3,  concentrates  on  the  eolutlon  of 
mixing  and  reactions  for  s  given  flowfleld. 
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Ttenfor*,  buoyancy  of  facts  end  the  Influence  of 
chemical  reaction  boat  ralaaae  on  the  flowfleld 
are  not  considered.  However ,  the  technique  avoids 
a  one  of  the  w  srlcal  problems  encountered  In  the 
full  numerical  a luula t lone.  The  comparisons  of 
the  conputed  probability  density  functions  with 
the  experimentally  determined  pdfs  show  some 
encouraging  agreements  for  the  reacting  case.  For 
the  non- reacting  case,  the  computed  pdf  Is 
antisymmetric  whereas  the  measured  pdf  shows  a 
mixing  bias  of  the  high-speed  stream  of  the  mixing 
layer.  A  complete  parameter  survey  of  Schmidt 
number,  Reynolds  number,  and  time  Is  underway. 


The  research  reported  In  this  paper  was  sponsred 
by  Air  Force  Office  of  Scientific  Research 
(AFOSR) . 
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Figure  1.  Infinite  Vortex  Row  Configuration. 


Figure  2.  Distortion  of  sn  Initially  horizontal 
Interface  Re  ■  50,  t  ■  0.02. 
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Figure  11.  Probability  density  of  reacting  0, 
species  for  Re  ■  125,  r  ■  0.004  and 
Sc  -  1.0. 
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Figure  12.  Probability  density  of  conserved 

scalar  for  Re  «  125,  t  ■  0.004  and 
Sc  -  1.0.  ^ _ . 


Figure  13.  Probability  density  for  the  non-react¬ 
ing  shear  layer  of  Masutanl  &  Bowaan 
at  x  -  lies.  0,  la  In  the  high  speed 
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Figure  14.  NO  aass  fraction  tor  tbe  simulated  shear 
layer  of  Masutanl  and  Bowman.  NO  is  in 
the  high  speed  stream  (upper  half-plane) 
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Figure  15.  Probability  density  of  0^  measured  in 
the  reacting  shear  layer  of  Masutanl 
and  Bowman  for  a  molar  0.  to  NO 
stoichiometry  of  .118.  "’no  Is  In  the 
high-speed  stream. 


Figure  16. 
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Probability  density  of  0.  In  the 
simulated  Masutanl  6  Bowman  experiment 
(conditions  of  Figure  15). 
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